Impurity Quantum Phase Transitions 
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We review recent work on continuous quantum phase transitions in impurity models, both with 
fermionic and bosonic baths - these transitions are interesting realizations of boundary critical 
phenomena at zero temperature. The models with fermion bath are generalizations of the standard 
Kondo model, with the common feature that Kondo screening of the localized spin can be suppressed 
due to competing processes. The models with boson bath are related to the spin-boson model of 
dissipative two-level systems, where the interplay between tunneling and friction results in multiple 
phases. The competition inherent to all models can generate unstable fixed points associated with 
quantum phase transitions, where the impurity properties undergo qualitative changes. Interestingly, 
certain impurity transitions feature both lower-critical and upper-critical "dimensions" and allow for 
epsilon-type expansions. We present results for a number of observables, obtained by both analytical 
and numerical renormalization group techniques, and make connections to experiments. 



I. INTRODUCTION 

Quantum mechanical systems can undergo zero- 
temperature phase transitions upon variation of a non- 
thermal control parameteriS, where order is destroyed 
solely by quantum fluctuations. Quantum phase transi- 
tions occur as a result of competing ground state phases, 
and can be classified into first-order and continuous tran- 
sitions. The transition point of a continuous quantum 
phase transition, the so-called quantum-critical point, is 
typically characterized by a critical continuum of excita- 
tions, and can lead to unconventional behavior - such as 
non-trivial power laws or non-Fermi liquid physics - over 
a wide range of the phase diagram. 

An interesting class of quantum phase transitions are 
so-called boundary transitions where only the degrees of 
freedom of a subsystem become critical. In this paper 
we consider impurity transitions - the impurity can be 
understood as a zero-dimensional boundary - where the 
impurity contribution to the free energy becomes singular 
at the quantum critical point. Historically, the first tran- 
sitions considered were the ones in the dissipative spin- 
bosonSiii and the anisotropic Kondo models^'^. Impurity 
quantum phase transitions require the thermodynamic 
limit in the bath system, but are completely independent 
of possible bulk phase transitions in the bath. 

In this paper we review recent work on quantum im- 
purity models undergoing zero-temperature phase transi- 
tions, associated with intermediate-coupling fixed points. 
These transitions can occur both in Kondo-type models 
with a fermionic bath and in quantum dissipative models 
with a bosonic bath. We shall show that these two classes 
are not fundamentally different, and the transitions - al- 
beit in different universality classes - can be analyzed us- 
ing similar techniques. We elaborate on the existence of 
both lower-critical and upper-critical "dimensions" and 
the possibility to access critical behavior via epsilon-type 
expansions. 

All models have the general form 



(1) 



Hh contains the bulk degrees of freedom, which gener- 
ically are interacting; however, under certain circum- 
stances the self-interaction is irrelevant and can be dis- 
carded from the outset. Tiimp contains the impurity de- 
grees of freedom, e.g., one or more quantum spins, to- 
gether with their coupling to the bath, which typically is 
local in space. Non-trivial quantum critical behavior in a 
model of the form obtains only if the thermodynamic 
limit for the bulk system Ti.^ is taken before the T ^ 
limit. 

The paper is organized as follows: In Sec. ^ we out- 
line possible types of transitions which have distinct 
finite-temperature crossover behavior, namely first-order, 
second-order, and infinite-order transitions. Issues re- 
lated to the bath Hh are discussed in Sec. IIIII in par- 
ticular we emphasize that in certain cases bosonic baths 
cannot be treated as Gaussian. Sec. IIVI introduces im- 
portant observables, together with their scaling behavior 
and related critical exponents. The remainder of the pa- 
per deals with concrete model systems showing impurity 
transitions. We start with the most prominent example 
of intermediate-coupling impurity physics, found in the 
fermionic multi-channel Kondo problem (Sec. 01. Dis- 
tinct non-trivial behavior can arise from a power-law vari- 
ation of the bath density of states at low energies - this 
is described for the pseudogap Kondo model (Sec. IVI|I . 
the sub-ohmic spin-boson model (Sec lVIll) . and the Bose 
Kondo model (Sec lVllTll . Impurities with both fermionic 
and bosonic baths are subject of Sec. lIXI whereas Secs.lXI 
and lXII deal with systems of coupled and multi-orbital im- 
purities, respectively. Finally, Sec. IXIlj is devoted to the 
quantum-classical mapping commonly employed to dis- 
cuss quantum phase transitions. We highlight that this 
mapping is inapplicable for many situations encountered 
in the context of quantum impurity models. A discus- 
sion of general aspects will close the paper. Most results 
quoted in this review have been obtained by large-A'^, 
perturbative renormalization group (RG) and Wilson's 
numerical renormalization group'^'^ (NRG) methods; we 
shall also mention some results from Bethe ansatz and 
conformal field theory techniques. Applications from var- 
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FIG. 1: Schematic finite-temperature phase diagrams for impurity quantum phase transitions; A and B are the stable phases 
(e.g. the localized and delocalized phases in the spin-boson model), a) First-order transition, i.e., level crossing. The finite- 
temperature properties are simply a thermodynamic mixture of A and B. b) Second-order transition, with power laws. The 
quantum-critical region, bounded by T* oc \t\'^ (where t = {r ~ r^jvc is the dimensionless measure of the distance to the tran- 
sition point), is controlled by an unstable RG fixed point, c) Infinite-order transition of Kosterlitz-Thouless type. No unstable 
fixed point is present, and the leading thermodynamic behavior displays only one crossover line. T* vanishes exponentially at 
the transition point. 



ious fields, e.g., magnetic impurities in non-metallic sys- mesoscopic systems, will be mentioned throughout the 
tems, strongly interacting quantum dots, and noise in text. 

I 



II. TRANSITIONS, CRITICAL DIMENSIONS, 
AND FINITE-TEMPERATURE CROSSOVERS 

Impurity quantum phase transitions can be continu- 
ous or discontinuous (first-order); the class of continu- 
ous transitions can be subdivided into transitions with 
power-law behavior (typically second-order) and those 
with exponential behavior (infinite-order transition, e.g., 
Kosterlitz-Thouless). Quantum impurity models permit 
sharp boundary transitions only at T = 0. Depending on 
the type of transition, the finite-temperature crossovers 
are different, see Fig.^ 

A first-order transition is a simple level crossing in the 
system's ground state: There exist two disconnected min- 
ima in the energy landscape of the system, and the finitc- 
T properties above the transition point are a simple ther- 
modynamic mixture of the two phases. A true quantum 
critical region is present for second-order transitions - 
this will be the main focus of the present article. 

Some of the examples presented below show second- 
order transitions with exponents depending upon a con- 
tinuous parameter, which specifies the exponent of a low- 
energy power law in the bath density of states which 
takes the role of a "dimension" - this applies to the 
pseudogap Kondo, the spin-boson, and the Bose-Fermi 
Kondo problems. As a function of this "dimension", 
the line of second-order transitions can terminate at 
a "lower-critical dimension" , with diverging correlation 
length exponent as this dimension is approached. Pre- 
cisely at the lower-critical dimension the transition can 



turn into one of Kosterlitz-Thouless type. The pseudogap 
Kondo problem also features an "upper-critical dimen- 
sion" where the interacting critical fixed point turns into 
a non-interacting one. More precisely, above the upper- 
critical dimension the transition becomes first-order with 
perturbative corrections, see Sec. 



III. BATHS 

The baths appearing in quantum impurity models typ- 
ically consist of fermions or bosons; in the latter case, 
these can be either real particles or quasiparticle collec- 
tive excitations. Of relevance for the impurity behavior 
is the local density of states of the bath at the impu- 
rity location, with an ultraviolet cutoff A and a certain 
asymptotic low-energy form. Most interesting here are 
gapless spectra with power-law behavior at low energies, 
as will appear in many models below. 

However, in a number of cases properties of the bath 
beyond its local spectrum are important, and this must 
be discussed in the context of the concrete physical model 
underlying the bath. In particular, a free Bose field with 
a gapless spectrum is a delicate object which can become 
unstable under infinitesimal perturbations (this should 
be contrasted from a free Fermi field, which has a ro- 
bust stability). One instance of such an instability is the 
response of spin-1 bosons to an applied magnetic field 
where the system is unstable to arbitrarily large fiuctu- 
ations. This means that the self-interaction of the bath 
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particles cannot be neglected, and must be treated at an 
equal footing with coupling between impurity and bath. 
Technically, there is a non-trivial "interference" between 
the two couplings, and the bath self- interaction signifi- 
cantly modifies the environment coupling to the impu- 
rity. It is thus not permissible to treat the environment 
as a Gaussian quantum noise, and focus exclusively on 
the coupling to the impurity. This applies e.g. to mag- 
netic fluctuations at a bulk quantum critical point below 
its upper-critical dimension, as exemplified for the Bose 
Kondo model in Sec. IVIIII 

Conversely, in models with non-interacting bath par- 
ticles the environment is Gaussian, and progress can be 
made by integrating out the bath. This results in a quan- 
tum model for the impurity only, with a (typically long- 
range) self-interaction in imaginary time direction. Such 
a model is effectively (0-1- l)-dimensional, and is believed 
to be asymptotically equivalent to a one-dimensional 
classical (spin) model by the virtue of the quantum- 
classical mapping. A paradigmatic example is the l/r^ 
Ising chain which is equivalent to both the ohmic spin- 
boson model and the anisotropic Kondo modcl'^. How- 
ever, recent investigations of the sub-ohmic spin-boson 
model indicate that the naive quantum-classical mapping 
can fail for certain types of long-range interactions^, see 
Sec. IXIII For spin impurity models with higher symme- 
try [XY or SU(2)] the quantum-classical mapping also 
fails; this can be traced back to the effect of the Berry 
phase describing the spin dynamics, which has no classi- 
cal analogue. 

IV. OBSERVABLES AND SCALING 

In this section we introduce a few observables which 
can be used to characterize the phases and phase tran- 
sitions discussed below. The equations will be given for 
a situation where the impurity degree of freedom is a 
single SU(2) spin S of size S; generalizations to other 
impurities are straightforward. The models to be dis- 
cussed display both "screened" and "unscreened" phases, 
where a "screened" phase has a unique ground state with 
quenched impurity degrees of freedom, whereas an "un- 
screened" phase shows a degenerate ground state. 

A. Susceptibilities 

Canonical response functions can be derived using a lo- 
cal magnetic field coupled to the impurity spin, by adding 

-^imp,aSa (2) 

to the impurity Hamiltonian Tiimp- For a spinful bath, an 
external magnetic field can also be applied to the degrees 
of freedom in Tib- For conduction electrons Ca-{x) the 
appropriate form is 

-HUx){ci<J^.,c^'){x) (3) 



The bulk field varies slowly as function of the space 
coordinate. 

With these definitions, a spatially uniform field applied 
to the whole system corresponds to — i?imp — H. 
Response functions can be defined from second deriva- 
tives of the thermodynamic potential, F — T In Z, in 
the standard wayi^^: u rneasures the bulk response 
to a field applied to the bulk, Xhnp imp the impurity 
response to a field applied to the impurity, and Xu,?mp 
is the cross-response of the bulk to an impurity field. In 
the absence of global SU(2) symmetry, response functions 
for different field directions have to be distinguished, and 
only some of them will be singular at criticality. 

The impurity contribution to the total susceptibility is 
defined as 

Ximp(r) 

— Ximp.imp "I" 2Xu,imp + {x u,u Xu.u 1, (4) 

where x^i^u^ is the susceptibility of the bulk system in 
absence of the impurity. For an unscreened impurity spin 
of size S we expect Ximp(T ^ 0) = S{S + 1)/(3T) in 
the low-temperature limit - note that this unscreened 
moment will be spatially smeared out due to the residual 
coupling to the bath. A fully screened moment will be 
characterized by Tximp ~ 0. In the presence of global 
SU(2) symmetry, the susceptibility Ximp does not acquire 
an anomalous dimension at criticalityi^i (in contrast to 
Xioc below), because it is a response function associated 
to the conserved quantity 5tot- Thus we expect a Curie 
law 

Cim 

^hn Ximp(r) = , (5) 

where the prefactor Cimp is in general a non-trivial uni- 
versal constant different from the free-impurity value 
S'(S' -I- 1)/3. Apparently, Eq. 10) can be interpreted as the 
Curie response of a fractional effective spinifl - examples 
are e.g. found in the pseudogap Kondo model (Sec. IVI|I 
and in the Bose Kondo model (Sec. IVIII|I below. 
The local impurity susceptibility is given by 

Xloc 

(T) 

— Ximp, imp j (6) 

which is equivalent to the zero-frequency impurity spin 
autocorrelation function. In an unscreened phase we have 
Xloc oc 1/T as r — !■ 0. This Curie law defines a residual 
local moment toioc at T = 0, which is the fraction of 
the total, freely fluctuating, moment of size S, which 
remained localized at the impurity site: 

2 

^unxioc(T) = ^. (7) 

A decoupled impurity has m^^^ — Cimp — S{S + l)/3, 
but a finite coupling to the bath implies mf^^ < Cimp. 
The quantity mioc turns out to be a suitable order 
parameter— for the phase transitions between an un- 
screened and a screened spin: at a second-order transition 
it vanishes continuously as f — > 0~. Here, t = (r — rc)/rc 
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is the dimensionless measure of the distance to criticahty 
in terms of couphng constants, with t > {t < 0) plac- 
ing the system into the (un)screened phase. Thus, Txioc 
is not pinned to the value of S{S + l)/3 for t < (in 
contrast to Tximp)- 

For models where the bulk does not carry spin, like 
the spin-boson model, a local susceptibility can still be 
defined, and X\oc ^'^^^ be singular at criticality. 



B. Scaling ansatz and critical exponents 



which describes critical local-moment fluctuations, and 
the local static susceptibility follows 

Xioc(T) = ^$2(^^j . (13) 

which reduces to the form quoted in (|10|) at criticality, 
t = 0. Here, 77-^ = 1 — x is a universal anomalous ex- 
ponent, and $1^2 are universal crossover functions (for 
the specific critical fixed point), whereas Bi^2 are non- 
universal prefactors. 



A scaling ansatz for the impurity part of the free energy 
takes the form 



(8) 



where t measures the distance to criticality, and iJimp 
is the local field, v is the correlation length exponent 
which describes the vanishing of the energy scale T*, 
above which critical behavior is observed^: 



T* oc \t\ 



(9) 



Note that there is no independent dynamical exponent z 
for the present (O-l-l)-dimensional models, formally z — \. 
The ansatz ^ assumes the fixed point to be interacting; 
for a Gaussian fixed point the scaling function will also 
depend upon dangerously irrelevant variables. 

With the local magnetization Mioc = {Sz) ~ 
-dFimp/dHitnp - note that Afioc(T ^ 0) = mioc 
O - and the corresponding susceptibility xioc ~ 
— 5^Fimp/(977iinp)^ we can define critical exponents as 
usual: 



(10) 



Mioc(t <0,T = 0,i/i„,p^O) cx i-tf, 

Xioc(t>0,r = 0) cx t-\ 

Afioc(t = 0,r = 0) cx \Hi^p\^/^, 

Xioc(i = 0,r) K T-^ 

xL(i = 0,r = 0,c.) cx |a;r^sgn(c.). 



The last equation describes the dynamical scaling of the 
local susceptibility. 

In the absence of a dangerously irrelevant variable, 
there are only two independent exponents. The scaling 
form (jHl allows to derive hyperscaling relations: 



C. Impurity entropy 

In general, zero-temperature impurity critical points 
can show a non-trivial residual entropy [contrary to bulk 
quantum critical points where the entropy usually van- 
ishes with a power law, S{T) oc T^]. The impurity en- 
tropy is evaluated as entropy of the system with impurity 
minus entropy of the bath alone. The stable phases usu- 
ally have impurity entropies of the form S'imp(r — > 0) = 
lug where is the integer ground state degeneracy, e.g., 
g = 1 for a Kondo-screened impurity and g = 2S + 1 for 
an unscreened spin of size S". At a second-order transition 
g can take fractional values; here it is again important 
that the thermodynamic limit in the bath is taken before 
the T hmit. 

Thermodynamic stability requires that the total en- 
tropy of a system decreases upon decreasing tempera- 
ture, dTS{T) > 0. This raises the question of whether 
the impurity part of the entropy, S'imp, has to decrease 
under RG flow (which is equivalent to decreasing T). 
The so-called g-theoremi^ provides a proof of this conjec- 
ture for systems with short-ranged interactions; for most 
quantum impurity problems this appears to apply. In- 
terestingly, both the pseudogap Kondo model and the 
spin-boson model provide explicit counter-examples , see 
Sees. IVII and IVIII with details in Refs 
another counter-example see Ref. T7.) 

On the technical side, we note that the impurity part of 
the thermodynamic potential, Fi,„p, will usually diverge 
with the ultraviolet cutoff A of the bath. Thus we have 
-Fimp — Emp — TSimp, where Simp is the non- universal 
(cutoff-dependent) impurity contribution to the ground- 
state energy. However, the impurity entropy S'imp is fully 
universal in the limit T/A — > 0. 
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(3 = j——, 2/3 + 7^ J/, 5 



2x 



1-x 



(11) 



D. 



T matrix 



Furthermore, hyperscaling also implies x = y. This is 
equivalent to so-called uj/T scaling in the dynamical be- 
havior - for instance, the local dynamic susceptibility will 
obey the full scaling formi 



OJ 



$1 - 



rpl/v 



T' t 



(12) 



Depending on the type of the impurity, further ob- 
servables can be used to characterize the behavior near 
criticality. In an Anderson or Kondo model the conduc- 
tion electron T matrix describes the scattering of the 
conduction electrons off the impurity. For an Anderson 
model, the T matrix is (up to a prefactor) equivalent to 
the full impurity electron Green's function. As with the 
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local susceptibility, we expect a scaling form of the T ma- 
trix spectral density near the intermediate-coupling fixed 
points similar to Eq. H12|) . In particular, at criticality a 
power law occurs: 



T{ljj) oc 



1 



,1-J7t 



(14) 



We will evaluate the anomalous exponent rjT for certain 
models below. 

Notably, the T matrix can be directly observed in ex- 
periments on Kondo impurities, due to recent advances 
in low-temperature scanning tunneling microscopy. This 
has been demonstrated for magnetic ions on metal 
surfacesiSiiS, and also with impurity moments in high- 
temperature superconductorsSfl.. Further, both metallic 
and semiconductor quantum dots in the Coulomb block- 
ade regime can show Kondo physics^^^, and transport 
experiments through the dot then probe the local spec- 
tral function via the differential conductance. 



E. Local non-Fermi liquid behavior 

Magnetic impurities in metals which show the standard 
Kondo effect can be described in terms of local Fermi 
liquid physics, i.e., the imaginary part of the conduction 
electron self-energy follows uP' at low energies, the impu- 
rity contribution to the entropy is etc. In general, the 
impurity influence on the bulk system can be captured by 
a potential scatterer with energy-dependent phase shift. 

In contrast, the intermediate-coupling fixed points to 
be described below are associated with local non-Fermi 
liquid behavior, and the above concepts do not apply. 
This is immediately clear e.g. for the impurity entropy 
which can be finite as T — *■ 0. A paradigmatic and well- 
studied example for local non-Fermi liquid behavior is 
found in the two-channel Kondo effect fSec. IVT). 



V. MULTI-CHANNEL KONDO MODEL 

Many of the models to be discussed in this paper are 
built from magnetic moments which can show the Kondo 
effect®. Originally, this effect describes the behavior of 
localized magnetic impurities in metals. The relevant 
microscopic models are the Kondo model and the single- 
impurity Anderson model. In the standard case (i.e. a 
single magnetic impurity with spin 1/2 coupled to a sin- 
gle conduction band with a finite density of states (DOS) 
near the Fermi level) screening of the magnetic moment 
occurs below a temperature scale Tk- The screening is 
associated with the flow to strong coupling of the effec- 
tive interaction between impurity and host fermions, i.e., 
perturbation theory diverges on the scale Tk- The Kondo 
temperature Tk depends exponentially on the model pa- 
rameters. For T <C Tk the system behaves as a local 
Fermi liquid. 



Kondo screening is strongly modified if two or more 
fermionic screening channels compete and the moment is 
then overcompensated. Nozieres and BlandinS^ proposed 
a two-channel generalization of the Kondo model, which 
shows overscreening associated with an intermediate- 
coupling fixed point and non-Fermi liquid behavior in 
various thermodynamic and transport properties. 

The Hamiltonian for a multi-channel Kondo model can 
be written as H = Hmk + X^i i' '^i^h 



K 



Hmk - ^J.S-s,(0) 



(15) 



H 



(f) 



-A 



where i — I, . . . , K is the channel index, and summation 

over repeated spin indices a is implied. The fermionic 

(f) 

baths Tijj are represented by one-dimensional chiral 
fermions Cfc, with a density of states po at the Fermi 
level and an ultraviolet (UV) cutoff A. The impurity 
spin S is coupled to the conduction electron spins at site 

0, Sa,(0) = 4(0X,,C,,,(0)/2 with C,,(0) = / dkkCka^, 

and cr" is the vector of Pauli matrices. 

In general, overscreening occurs for any number of 
channels K > 1 coupled to a spin 1/2. It does not require 
fine-tuning of the Kondo coupling, provided that the cou- 
plings Ji in different channels are equal, i.e., if the system 
obeys SU(2)spinX SU(ii:)channei symmmctry^. A quan- 
tum phase transition can be driven by a channel asym- 
metry, with the RG fiow diagram for the two-channel 
case shown in Fig. [21 One of the two equivalent strong- 
coupling Fermi liquid fixed points is reached for un- 
equal couplings Ji ^ J2, whereas for Ji = J2 the sys- 
tem fiows to the two-channel intermediate-coupling fixed 
point. Thus, the two-channel fixed point acts as critical 
fixed point separating the two Fermi liquid phases. 

As an aside we note that the two-channel Ander- 
son model has recently been shown to have even richer 
behavior—: it displays a line of non-Fermi liquid fixed 
points which continuously connects the two-channel spin 
and two-channel charge Kondo effects. 



RG 



The existence of an intermediate-coupling fixed point 
for a large number of screening channels K can be seen 
from the RG fiow equation for the dimensionless Kondo 
coupling, j = ji = j2 {ji = poJi initially), in the channel- 
symmetric case 



(16) 



to two-loop order. It predicts an unstable fixed point 
at j* — 2/K, which indeed exists for all K > 2 and 
is perturbatively accessible in the limit of large channel 
number, K > ]Mi^. 
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LM 



We caution the reader that the quoted behavior of 
Ximp violates the naive scahng expectation Ximp l/T 
(0) . Even though the low-energy physics is controlled by 
an intermediate-coupling fixed point, a "compensation" 
effectil causes most thermodynamic response functions 
to vanish in the naive scaling limit, and the leading low- 
temperature behavior is exposed only upon considering 
corrections to scahng^^-^-. For this reason Ximp and xioc 
follow the same anomalous behavior, in contrast to the 
pseudogap Kondo and Bose Kondo models, described in 
Sees. El and Eml 



D. Applications 



FIG. 2: RG flow diagram for the two-channel Kondo model 
where the two axes denote the renormalized couplings ji,j2 
in the two channels, with ji = poJi initially. LM is the 
local-moment fixed point of a decoupled impurity, ICK is the 
strong-coupling fixed point of the single-channel Kondo ef- 
fect, and 2CK is the intermediate-coupling two-channel fixed 
point. 

B. Exact solution 

The multi-channel Kondo model can be exactly solved 
by means of the Bethe ansataSl. In addition, many of 
the low-energy properties of the two-channel Kondo and 
related models have been successfully studied using con- 
formal field theory techniques'^. This allows to calculate 
the leading temperature dependencies of physical observ- 
ables, like entropy, suceptibility, and also dynamic mul- 
tipoint correlation functions which are not directly ac- 
cessible to Bethe ansatz methods. For a more detailed 
description we refer the reader to the Iiterature22i22i2£. 

C. Observables 

The anomalous low-temperature properties at the 
multi-channel non-Fermi liquid fixed point are controlled 
by the leading irrelevant coupling (see below), which has 
scaling dimension y = 2/(2-1- K). One finds for the im- 
purity susceptibility and the specific heat ratio 

Ximp « {T/T^fy-\ 

l = Q^p/T <x [T/TKfy-K (17) 

In the two-channel case, K — 2, logarithmic corrections 
occur, leading to 

Ximp oc \n{TK/T) , 
7 = amp/T (X HTk/T). (18) 

For K = 2 the residual entropy is S'imp — 5 In 2, and 
an anomalous Wilson ratio R = x/l = 8/3 appears, 
in contrast to the result for the standard Kondo model 
i? = 8/4 = 2. 



Experimental realizations of two-channel Kondo 
physics have been discussed early on in the context of 
certain rare-earth compounds. The idea is that in ma- 
terials like UBei3 orbital (i.e. quadrupolar) degrees of 
freedom can be screened by the interaction with conduc- 
tion electrons, and the electron spin a =f , i provides the 
two channels required for overscreening. This proposal 
has been controversial, and we refer the reader to a re- 
cent review24 for details. 

As the Kondo effect of a single spin can be conveniently 
investigated in transport experiments through quantum 
dots in the Coulomb blockade regim e- - i ' ^^ , various pro- 
posals for the realization of two-channel Kondo physics 
in quantum-dot devices have been put forward. The dif- 
ficulties for the spin two-channel Kondo effect are that 
the screening channels need to be independent (i.e. no 
mixing due to cotunneling between the bands) and the 
couplings to the channels need to be equal. Provided that 
cotunneling is suppressed, the channel asymmetry may 
be tuned by suitable gate voltages, leading to a quantum 
phase transition as illustrated in Fig. |21 the transport 
signatures of this phase transition have been detailed in 
Ref. lU Versions of the two-channel Kondo effect utiliz- 
ing charge degrees of freedom in quantum dots'^^ natu- 
rally satisfy the above requirements, as the two indepen- 
dent channels are provided by the electron spin. Here, 
the problem is to reach temperatures sufficiently below 
Tk- 

A recent proposal by Oreg and Goldhaber-GordonSi 
utilizes Coulomb blockade physics in the leads, formed 
by metallic islands, to suppress cotunneling between dif- 
ferent leads. Subsequently it has been shown—- that a 
tiny charging energy in those leads is sufficient to sta- 
bilize multi-channel physics as T — > 0. In such a situ- 
ation, a single-channel Kondo resonance develops at in- 
termediate temperature, which is destroyed at lowest T, 
i.e., the system crosses over from single- to multi-channel 
physics (!) upon lowering temperature. This "reversed" 
renormalization group flow is opposite to the standard 
behavior in a multi-channel model with small channel 
anisotropy. Fig. |21 It contradicts the conventional wis- 
dom that single-channel Kondo physics cannot be de- 
stroyed by small perturbations - the physical reason is 
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the long-range nature of the Coulomb charging energj*ii. 

We note that to date experimental efforts in realizing 
the two-channel Kondo effect in quantum dots have re- 
mained unsuccessful. 

VI. PSEUDOGAP KONDO AND ANDERSON 
MODELS 

A straightforward possibility to suppress Kondo 
screening in a single-channel situation is to reduce the 
electron bath density of states at the Fermi level to zero. 
The absence of low-energy states prevents screening for 
small Kondo couplings. If the bath gap arises from host 
superconductivity this can be interpreted as competition 
between Kondo singlet formation and Cooper pairing. 

Two different types of gaps have to be distinguished: 
so-called hard-gap and pseudogap (soft-gap) systems. In 
the hard-gap case, the DOS p{e) is zero in a finite energy 
interval around the Fermi level - this is e.g. realized in 
a conventional superconductor. The resulting impurity 
transition between a local-moment (LM) phase without 
Kondo screening, realized at small Kondo coupling J, 
and a screened strong-coupling (SC) phase, reached for 
large J, is of first order, and it occurs only in the presence 
of particle- hole (p-h) asymmetry25i. In the p-h symmetric 
case, the local-moment state persists for arbitrary values 
of the coupling. 

The pseudogap case, first considered by Withoff and 
Fradkin'^®, corresponds to a bath with p{e) oc \e\^ (r > 0), 
i.e., the DOS is zero only at the Fermi level. The cor- 
responding Kondo and single-impurity Anderson models 
interpolate between the metallic case (r = 0) and the 
hard-gap case (r oo). The pseudogap case < r < oo 
leads to a very rich behaviour, in particular to a con- 
tinuous transition between a local-moment and a strong- 
coupling phase^Si^LSi^ Central to our discussion will be 
the single-impurity Anderson model with a pseudogap 
host density of states, H = Ha + 'Hh^' 

Uk = eo/i/. + C/ori/Tn/i+yo(/^c,(0) + h.c.)(19) 
= / dk\k\' kcl^cu. 

J -A 

with notations as in Sec. and the bath is now rep- 
resented by linearly dispersing fermions in (1 -I- r) di- 
mensions. In the Kondo limit of the Anderson model 
charge fluctuations are frozen out, and the impurity site 
is mainly singly occupied. Via Schrieffer- Wolff transfor- 
mation one obtains the standard Kondo model, H. — 
nK+nl^\ with 

Hk = JkS ■ s(0) . (20) 

The Kondo coupling is related to the parameters of the 
Anderson model H19|) through: 
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FIG. 3: T = phase diagram for the pseudogap Ander- 
son model in the p-h symmetric case (solid line, f/o = 10"'^, 
£o = —0.5 ■ 10~^, conduction band cutoff A = 1) and the p-h 
asymmetric case (dashed line, eo = —0.4 • 10~^); A mea- 
sures the hybridization strength, -rVq piio) = A|aj|''. The 
lines of critical points separate the strong coupling phase SC 
(A > Ac) from the local-moment phase LM (A < Ac). After 
Ref.E^ 

The Kondo limit is reached by taking J/q ~* oo? ^ 
— cx), Vo — > c», keeping Jk fixed. In the absence of p-h 
symmetry the Schrieffer- Wolff transformation also gen- 
erates a potential scattering term in the effective Kondo 
model^. 



A. Pseudogap Anderson model: Phase diagram 
and NRG 

Fig. O shows a typical phase diagram for the pseu- 
dogap Anderson model22ii£. In the p-h symmetric case 
(solid) the critical coupling A, measuring the hybridiza- 
tion between band electrons and local moment, diverges 
at r — 1/2, and no screening occurs for r > 1/2. No 
divergence occurs for p-h asymmetry (dashed). Impor- 
tantly, the strong-coupling phases in the p-h symmetric 
and asymmetric situations belong to different low-energy 
fixed points, in the following denoted by SSC and ASC, 
respectively. 

We now briefly describe the stable phases; NRG cal- 
culations have established that their low-energy prop- 
erties are identical for the pseudogap Anderson and 
Kondo models'^^ . Due to the power-law conduction band 
DOS, already the stable LM and SSC/ASC fixed points 
show unconventional behaviori224£. The LM phase has 
the properties of a free spin 1/2 with residual entropy 
Simp = In 2 and low-temperature impurity susceptibil- 
ity Ximp = l/(4r), but the leading corrections show r- 
dependent power laws. The p-h symmetric SSC fixed 
point has very unusual properties, namely S'imp = 2r In 2, 
Ximp = r/{8T) for < r < 1/2. In contrast, the p- 
h asymmetric ASC fixed point simply displays a com- 
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pletely screened moment, Simp = T'ximp = 0. The impu- 
rity spectral function follows a lo^ power law at both the 
LM and the ASC fixed point, whereas it diverges as lj^^ 
at the SSC fixed point - this "peak" can be viewed as 
a generalization of the Kondo resonance in the standard 
case (r — 0), and scaling of this peak is observed upon 
approaching the SSC-LM phase boundarji^ML. 

Let us turn to the critical fixed points. Early work em- 
ployed a weak-coupling RG based on an expansion in the 
Kondo coupling, see next subsection. It turns out that 
the applicability of this approach is restricted to small r. 
Interestingly, the NRG studies^ showed that the fixed- 
point structure changes at r = r* « 0.375 and also at 
r = i , rendering the interesting case of r = 1 inaccessible 
from weak coupling. For r < r* the p-h symmetric and 
asymmetric models reach the same critical fixed point 
(SCR) which is p-h symmetric. In contrast, for r > r* a 
p-h asymmetric critical fixed point (ACR) exists, which 
then controls the phase transition in the p-h asymmetric 
model. The SCR fixed point ceases to exist for r > 1/2. 
In addition, the critical fluctuations in the p-h asymmet- 
ric case change their character at r — 1: whereas for r < 1 
the exponents take non-trivial r-dependent values and 
obey hyperscaling, exponents are trivial for r > 1 and 
hyperscaling is violatedi^^. 

Analytical understanding came with RG expansions 
applied to the pseudogap Anderson modeli^ which al- 
low to formulate the universal critical theories for r near 
1/2 and near 1. In the p-h symmetric case the line 
of non-trivial phase transitions terminates at two lower- 
critical dimensions (!), r = and r = 1/2. The following 
two subsections describe two RG expansions, with small 
parameters r and (1/2 — r) respectively, which access 
the same critical fixed point (SCR) and are expected to 
match. In the p-h asymmetric case an expansion can be 
done in the hybridization around the valence-fluctuation 
point of the Anderson model. Bare perturbation theory is 
sufficient for all r > 1; for r < 1 a perturbative RG pro- 
cedure is required to calculate critical properties, with 
the expansion being controlled in the small parameter 
(1 — r). In particular, this identifies r = 1 as the upper- 
critical dimension of the (asymmetric) pseudogap Kondo 
problem, and consequently observables acquire logarith- 
mic corrections for r — 1. The full RG flow diagrams in 
the variables of the Anderson model H19|) are displayed 
in Figs. 21 a-iid [3 for the p-h symmetric and asymmetric 
cases, respectively. The properties near criticality will be 
discussed in detail in Sec. IVflT below. 



r = 0. The flow of the renormalized Kondo coupling j is 
given by the beta function 

P{.l)^rj~f. (22) 

For r > there is a stable fixed point at j* = cor- 
responding to the local- moment phase (LM). An unsta- 
ble fixed point, controlling the transition to the strong- 
coupling phase, exists at j* = r, and the critical prop- 
erties can be determined in a double expansion in r and 
j. P-h asymmetry is irrelevant, i.e., a potential scatter- 
ing terms scales to zero according to (3{v) = rv, thus 
the above expansion captures the p-h symmetric critical 
fixed point (SCR). 

C. Pseudogap Anderson model: RG near r = 1/2 

For r near 1/2 the p-h symmetric critical fixed point 
moves to strong Kondo coupling, and the language of the 
p-h symmetric Anderson modeU^ becomes more appro- 
priate. Expanding in the renormalized hybridization v 
and on-site interaction u, the fiow equations 

1 — r 

P{v) = —v + v\ 

/3(u) = (l-2r).-fc|M^3 (23) 

capture various fixed points; notably the RG equation 
for V is exact to all ordersi^. The RG flow is depicted in 
Fig.H 

For r < 1/2 a stable fixed point is at v*"^ = (1 — r)/2, 
u* — 0, which represents nothing but the physics of 
the non-interacting resonant level model with a power- 
law density of states - interestingly, the impurity spin is 
not fully screened for r > 0, and the residual entropy 
is 2rln2. This fixed point is identical with the sym- 
metric strong-coupling fixed point (SSC) of Gonzalez- 
Buxton and Ingersent-22^; it becomes unstable for r > 1/2. 
Additionally, for r < 1/2 there is a pair of critical 
fixed points (SCR,SCR') located at w*^ = (1 - r)/2, 
M*2 = 7r2(l-2r)/[3(7r-21n4)]. These describe the tran- 
sition between an unscreened (spin or charge) moment 
phase and the symmetric strong-coupling (i.e. screened) 
phase. Importantly, both l|22|) and (|23|l capture the same 
critical fixed point. Using (|23|l . the expansion is per- 
formed around the strong-coupling fixed point (SSC), 
and is valid for u* <C 1, i.e., for 1/2 - r ^ 1. 



B. Pseudogap Kondo model: RG near r — 

For small values of the DOS exponent r, the phase 
transition in the pseudogap Kondo model can be accessed 
from the weak-coupling limit, using a generalization of 
Anderson's poor man's scaling. Power counting about 
the local moment fixed point (LM) shows that the Kondo 
coupling J is marginal at the lower-critical dimension 



D. Pseudogap Anderson model: RG near r — 1 

In the p-h asymmetric case a phase transition occurs 
for all r > 0; here the strong-coupling phase is maximally 
p-h asymmetric and corresponds to a fully screened spin 
for all r. Remarkably, for 0<r<r*?a0.375 p-h symme- 
try is dynamically restored at the phase transition, and 
the critical properties are described by the expansions 
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FIG. 4: Schematic RG flow diagrams for the particle-hole 
symmetric single-impurity Anderson model with a pseudogap 
DOS, p{u)) oc \ujy ■ The horizontal axis denotes the renormal- 
ized on-site level energy e (related to the on-site repulsion u by 
u = — 2e), the vertical axis is the renormalized hybridization 
V. The thick lines correspond to continuous boundary phase 
transitions; the full (open) circles are stable (unstable) fixed 
points, for details see text, a) r — 0, i.e., the familiar metallic 
case. For any finite v the fiow is towards the strong-coupling 
fixed point (SC), describing Kondo screening, b) 0<r<l/2: 
The local-moment fixed point (LM) is stable, and the transi- 
tion to symmetric strong coupling (SSC) is controlled by the 
SCR fixed point. For r 0, SCR approaches LM, and the 
critical behavior at SCR is accessible via an expansion in the 
Kondo coupling j. In contrast, for r 1/2, SCR approaches 
SSC, and the critical behavior can be accessed by expanding 
in the deviation from SCR, i.e., in £ = —u/2. c) l/2<r < 1: v 
is still relevant at u = 0. However, SSC is now unstable w.r.t. 
finite u. At finite v, the transition between the two stable 
fixed points LM and LM' is controlled by SSC (which is now 
a critical fixed point!), d) r>l: v is irrelevant, and the only 
transition is a level crossing (with perturbative corrections) 
occurring at v = m = 0, i.e., at the free-impurity fixed point 
(FImp). After Ref.IH 
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FIG. 5: Schematic RG flow diagrams for the maxi- 
mally particle-hole asymmetric pseudogap Anderson impu- 
rity model. The horizontal axis denotes the on-site impurity 
energy, s, the vertical axis is the fermionic coupling v, the 
bare on-site repulsion is fixed at = cxd. The symbols are 
as in Fig. |1] a) r* < r < 1: w is relevant, and the tran- 
sition is controlled by an interacting fixed point (ACR). As 
r — > r* « 0.375, p-h symmetry at the critical fixed point is dy- 
namically restored, and ACR merges into the SCR fixed point 
of Fig. 2]- this cannot be described using the RG of Sec. I VI Dl 
In the metallic r = situation, the flow from any point with 
V ^ is towards the screened singlet fixed point with e = oo. 
b) r > 1: V is irrelevant, and the transition is a level crossing 
with perturbative corrections, occuring at v = e = 0, i.e., the 
valence-fiuctuation fixed point (VFl). After Ref. Il5l 



discussed above. For r > r* there exists an additional 
critical fixed point (ACR) with finite p-h asymmetry, and 
its properties are best discussed using an infinite-C/o An- 
derson model— . Expanding around the e = 0, i; = 
limit, one finds the flow equations 

r,, \ I — r 3 , 
P{v) = —v+-v\ 

/5(e) = -e + Sv^e - ^2 , (24) 

with a critical fixed point w*^ = (1 — r)/3, e* = — (l-r)/3 
for r < 1, and w* = e* = for r > 1. The complete RG 
flow diagrams are in Fig. |S1 

Note that the flow is very similar to the one of a (f)^ 
model, with (1 — r) playing the role of (4 — d), as the An- 
derson model hybridization is marginal at r = 1 . Clearly, 
the critical fixed point is interacting for r < 1 (the ana- 
logue of the Wilson-Fisher fixed point), whereas for r > 1 
we have a level crossing with perturbative corrections 
(VFl, the analogue of the Gaussian fixed point). This 
justifies to identify r — 1 with the upper-critical dimen- 
sion of the problem. RG analysis has shown that all 
fixed points presented above are stable with respect to 
spin anisotropics. 
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E. Pseudogap Kondo model: large-A'^ 

A large-A'^ theory can be constructed by generalizing 
the spin symmetry of both quantum impurity and con- 
duction electrons from SU(2) to SU(-/V). An antisym- 
metric representation of the impurity spin uses auxiliary 
fermions fa {a — 1, ... , N): 

Sal3 = flflS - qSal3 (25) 

with a constraint fafa — 1^ enforced by a chemi- 
cal potential Aq; the value q — 1/2 corresponds to p-h 
symmetry and will be used in the following. 

The single-channel Kondo model takes the form 

Hk = -^/Ic.(0)4(0)/^ + XoflU ■ (26) 

In the limit of ^ cxd the action is dominated by a 
static saddle point where the field b conjugate to /qCq,(0) 
acquires an expectation value. We arrive at the non- 
interacting model 

Hk = - (6/Ic„(0) + h.c.) + XoflU + ^ (27) 
and the self-consistency equation 

^=#E(4(0)//5)- (28) 

Eqs. (|27I28|) constitute the well-known slave-boson mean- 
field approximation. Non-zero b signals Kondo screening. 
If 5 at r = then b will vanish continuously as 
a temperature Tk is approached from below - this de- 
fines the Kondo temperature. Note that the sharp finite- 
temperature transition is an artifact of the large- A'^ limit. 

For a metallic density of states, r = 0, the slave-boson 
method yields the correct exponential dependence of Tk 
on Jk. For r > a T = transition at a finite Jk is 
predicted. However, the critical properties of this transi- 
tion are not reproduced correctly; in particular at r = 1 
the slave-boson method yields an essential singularity of 
Tk near the transition point^LS (instead of the linear 
behavior, see above). 

An interesting alternative to the single-channel large- 
A^ limit is a multi-channel version, where the impurity 
is coupled to K screening channels2&. The Hamiltonian 
then reads 

Hmk = -^/Ic„,(0)4(0)/^ + XoPja (29) 

where i — 1,...,K is the channel index. Taking the 
N ^ oo limit with K — yields a dynamic saddle 
point, characterized by a time-dependent propagator of 
the bosonic field Bi conjugate to f^Ccyi{0). One finds the 
following integral equations for the fermionic and bosonic 
Green's functions G/(t) = -(T,/„(r)/t (0)}, Gb(t) = 

S/(t) = 74"(^)Gb(t), SB(r) = ^G^^H~T)Gf{T) , 

(30) 
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FIG. 6: Inverse correlation length exponent l/i/ obtained 
from NRG, at both the symmetric (squares) and asymmetric 
(triangles) critical points, together with the analytical RG 
results from the expansions in r (solid), in (i — r) (dashed), 
and in (1 - r) (dash-dot), eq. After Ref.[ll 



where Gg is the local Green's function of the fermionic 
bath, and the self-energies S/ and T,b are defined by: 

Gj^{iLu„) = iw„ + Ao - i;/(ia;„) (31) 
GsH^'^n) = ^-Ss(*J^„). (32) 

•JK 

In these expressions w„ — (2n + l)7r//? and Vn = 2m:/ f3 
denote fermionic and bosonic Matsubara frequencies. 
The third saddle point equation fixes the impurity "oc- 
cupation" and thus determines Aq: 

Gf{T = 0-) = qo. (33) 

Eqs. H30I31I32|I are identical in structure to the usual 
NCA equations^, but the physics is changed by the con- 
straint eq. (|33|l which keeps track of the choice of the 
impurity spin representation. 

The analysis of the NCA equations for the pseudogap 
case has been presented in Ref. As expected for a 
multi-channel Kondo model the exactly screened phase 
at large Jk is replaced by an overscreened one. For r > 
this overscreened phase is bounded by second-order tran- 
sitions to phases with either a free spin moment (LM, 
at small Jk) or a free channel moment (the latter only 
exists for p-h asymmetry and for large Jk). In both 
the overscreened phase and at criticality the auxiliary- 
particle propagators G/ and Gb follow power laws at 
low energy. The phase transitions display r-dependent 
exponents - naturally these differ from the ones of the 
SU(2) single-channel model. Various critical properties 
are reproduced correctly, e.g., the a;"'' divergence of the 
conduction electron T matrix at criticality, see below. In 
the p-h asymmetric case r = 1 is recovered as upper- 
critical dimension. 



11 




0.2 0.4 0.6 0.8 1 



r 

FIG. 7: Numerical data for the impurity susceptibility, 
TXimp, at both the symmetric (squares) and asymmetric (tri- 
angles) critical points, together with the renormalized per- 
turbation theory results from the expansions in r [solid], in 
(i - r) [dashed], and in (1 - r) [dash-dot]. After Ref.lm 



F. Results for observables 

Let us quote some important observables near criti- 
cality for the SU(2)-symmetric pseudogap Anderson and 
Kondo models - these results have been obtained using 
NRG and by the three e-type expansions together with 
renormalized perturbation theory. 

The correlation length exponent v jnj can be obtained 
from the perturbative RG by expanding the beta func- 
tions in the vicinity of the fixed point. The results are 

fr-^ + Olr^) SCR, r<l 

2-4r + C»([i -r]2) SCR, i-r<l . (34) 
r-f ©([1 -r]2) ACR, l-r< 1 

At the symmetric critical fixed point v diverges for both 
r 0+ and r 1/2^. For r > 1 the transition is a level 
crossing, formally v — \. Fig. |3 shows a comparison of 
these results for v with NRG results. Other exponents 
can also be evaluated perturbatively, see Ref . [l^ and the 
full set of exponents (|10() can be derived from hyperscal- 
ing relations which are valid for r < 1. 

In the quantum critical regime unconventional behav- 
ior corresponding to a fractional moment can be ob- 
served. As explained in Sec. IIVI the impurity entropy 
S'imp and the impurity Curie moment Tximp (Fig.Ej) are 
universal functions of r in the low-temperature limit^^'^. 

Let us briefly mention our result for the conduction 
electron T matrix, describing the scattering of the c elec- 
trons off the impurity. At criticality it follows a power law 
p4|l . and one finds the exact result r/T — 1 — r for r < 1, 
i.e., for all interacting fixed points considered above the 
T matrix follows T{uj) oc w"*". This perfectly agrees with 
NRG calculations'"^. Furthermore, at the upper-critical 
dimension we find ImT(w) cx 1/(lu \ logwj^) - this applies, 
e.g., to a Kondo impurity in a d-wave superconductor. 



G. Applications 

The Kondo effect in a non-metallic host is most con- 
veniently studied using superconductors. As above, both 
impurity spins embedded in a bulk material and inter- 
acting quantum dots are candidate realizations. A recent 
experimenljiS employed a carbon nanotube dot coupled 
to Nb leads, which can be switched between s-wave su- 
perconducting and normal states by a small magnetic 
field. A sharp crossover in the transport properties was 
found as a function of T/f/A, where Tk is the normal- 
state Kondo temperature of the dot and A the gap of the 
superconductor, consistent with the expected first-order 
transition within a hard-gap Kondo model. Quantum 
dots coupled to unconventional superconductors remain 
to be realized. Note that a complete description of the 
transport requires to account for multiple Andreev reflec- 
tions. Other suggestions for pseudogap Kondo physics 
include quantum dots coupled to either interacting one- 
dimensional electron liquidsi^ or to disordered meso- 
scopic conductors*^. 

Impurity moments embedded in bulk unconventional 
superconductors (or other systems with nodal order pa- 
rameters) represent another realization of pseudogap 
Kondo physics. Indeed, in c?-wave high- Tc superconduc- 
tors (r — 1 here) non-trivial Kondo-like behavior has 
been observed associated with the magnetic moments in- 
duced by non-magnetic Zn or Li impurities^. In under- 
doped superconducting samples the local susceptibility 
as measured by NMR on the impurity site displays a 
Curie law which turns into a Curie- Weiss law at optimal 
doping. This apparent transition from an unscreened to 
a screened moment is likely driven by a change of either 
the size of the superconducting gap, i.e., a change of the 
dimensionless coupling in the pseudogap Kondo model, 
or of increasing strength of the bulk spin fluctuations 
upon underdoping which tend to suppress Kondo screen- 
ing (see Sec. lIX|) . STM experiments^^ show a large peak 
at small bias in the local density of states near the impu- 
rity site - this has been interpreted as arising from the 
Kondo screening of the impurity- induced moment, and 
explicit calculations within a p-h asymmetric pseudogap 
Kondo model have been presented^S*^. 



VII. SPIN-BOSON MODEL 

After having described fermionic impurity models with 
intermediate-coupling fixed points, let us turn to systems 
of quantum impurities coupled to bosonic baths, which 
are an equally interesting model class. They have been 
first introduced in the context of the description of dis- 
sipative dynamics in quantum systems^. The simplest 
realization is the so-called spin-boson model, describing 
a two-level system coupled to a single bath of harmonic 
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oscillators: 



Ao 



3B 



(35) 

It describes a spin, tunneling between | f) and | |) via 
Aq, and being damped by the coupling to the oscillator 
bath. The coupling between spin and bath is completely 
specified by the bath spectral function 



J{uj) — vr ^ (w — Ui) . 



(36) 



Of particular interest are power-law spectra 

J{lo) = 2tt auol^^'Lo'' , < < t^c , s > -1 (37) 

where Wc is a cutoff, and the dimensionless parameter a 
characterizes the coupling or dissipation strength. 

The case of s = 1 corresponds to the well-studied ohmic 
spin-boson model^, which shows a Kosterlitz-Thouless 
quantum transition, separating a localized phase at a > 
ac from a delocalized phase at a < a^^. In the local- 
ized regime, the tunnel splitting between the two levels 
renormalizes to zero, i.e., the system gets trapped in one 
of the states | t) or | |), whereas the tunnel splitting stays 
finite in the delocalized phase. In the limit Aq ^ ujc the 
transition occurs at ac = 1. The super-ohmic situation, 
s > 1 is known to show weakly damped dynamics for any 
dissipation strength a, i.e., the two- level system is always 
delocalized in the zero-temperature limit, provided that 
Aq ^ 0. Much less is known about the sub-ohmic case 
(s < 1) - this turns out to be very interesting and will 
be discussed below. 



A. Spin-boson model: Phase diagram and NRG 

In the following we focus on the sub-ohmic exponent 
range, < s < 1. It has recently been established that a 
continuous quantum transition between a localized and 
a delocalized phase exists for all bath exponents < s < 
li&. We emphasize that this transition does not appear in 
the popular non-interacting blip (NIBA) approximation^. 

The two stable phases of the sub-ohmic spin-boson 
model are conventional: The delocalized phase has a 
unique ground state and a finite susceptibility x^^j i-e., 
the level splitting is finite in the low-temperature limit. 
In contrast, the localized phase has a residual entropy 
•Simp = ln2 and follows a Curie law of the form ||7|). 
For fixed s there exists a single critical fixed point with s- 
dependent exponents; as function of s these fixed points 
form a line which terminates in the Kosterlitz-Thouless 
transition point at s = l. The transition also disappears 
in the limit s — > 0^, and the model is always in a local- 
ized ground state for s < 0. Numerical results for the 
phase diagram are shown in Fig. |S1 

It is generally believed that the spin-boson model is 
equivalent to a one-dimensional classical Ising model 
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FIG. 8; Phase diagram for the transition between delocalized 
(a < Qc) and localized phases [a > ac) of the spin-boson 
model 1351 for bias e = and various values of Aq, deduced 
from NRG. After Ref.^S- 



with long-range interactions^iSiSi, specifically the bath 
with Lo'^ damping rate corresponds to an Ising interac- 
tion falling off as l/r^^". This quantum-classical map- 
ping will be discussed in more detail in Sec. IXIII while 
this equivalence holds for s > 1/2 including the ohmic 
case, it has recently been shownS that it is violated for 
s < 1/2. 



B. Spin-boson model: RG near s = 1 

We now turn to analytical approaches to describe the 
spin-boson phase transition. Power counting about the 
Aq = a — fixed point shows that the coupling a is 
marginal at s = 1, suggesting that a perturbative RG 
may be possible. However, for the spin-boson model, 
i.e., Ising symmetry, one encounters a flow of a towards 
strong coupling for s < 1. Progress can be made using 
a different formulation of the partition function, namely 
using a kink gas representation, where the kinks represent 
spin flips along the imaginary time axis or equivalently 
Ising domain walls---:. A perturbative RG analysis of this 
Ising model has been performed by KosterlitaS. This 
expansion, controlled by the smallness of the kink fugac- 
ity, is done around the ordered phase of the Ising model, 
corresponding to the localized fixed point of the spin- 
boson quantum problem. Carrying over these results to 
the spin-boson model, ones arrives at RG equations: 



/3(a) 
/3(A) 



a(A2- 
-A(l 



s-1), 
-a). 



(38) 



valid for small A, where A = Aq/oJc is the dimensionless 
tunneling strength. The RG flow is sketched in Fig. 1 of 
Rcf. 50. For s = 1, these equations are equivalent to the 
ones known from the anisotropic Kondo model, and de- 
scribe a Kosterlitz-Thouless transition, with a fixed line 
A = 0, a > 1. For s < 1 there is an unstable fixed 
point at a = 1, A^ = 1 — 5; clearly it is perturbatively 
accessible for small values of (1 — s) only. 
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FIG. 9; Numerical data^ for the correlation length exponent 
l/i/ obtained from NRG, together with the two perturbative 
RG results in Eq. (1401 which are asymptotically valid near 
s = 1 and s — 0, respectively. 



C. Spin-boson model: RG near s — 

A different RG expansion can be performed around 
the delocalized fixed point of tfie spin-boson model, cor- 
responding to finite Aq and infinitesimal a. For conve- 
nience we shall assume equal couplings, Xi = A, then the 
energy dependence of J{uj) is contained in the density of 
states of the oscillator modes Ui, and we have a cx A^. 
Eigenstates of the impurity are | —>x) and | <— a;), with 
an energy splitting of Aq. The low-energy Hilbert space 
contains the state | -^x) only, and interaction processes 
with the bath arise in second-order perturbation theory, 
proportional to kq — A^/Aq. Power counting shows that 
kq is marginal at s = 0. Performing RG within the 
low-energy sector one finds the RG beta function for the 
renormalized second-order tunneling amplitude: 

P{k) = sk-k^ . (39) 

Besides the stable delocalized fixed point k = this 
flow equation displays an infrared unstable fixed point 
at K* = s which controls the transition between the de- 
localized and localized phases. It is worth pointing out 
that no (dangerously) irrelevant variables are present in 
this theory. 

As an aside, we note that at s = the bath coupling 
is marginally relevant. Therefore the impurity is always 
localized as T — > 0, with a localization temperature given 
by T* =u>aeM-^ouJc/X'^)- 



D. Results for observables 

We now turn to results for various critical properties 
of the spin-boson model, determined using the two RG 
expansions and via bosonic NRG. 

The correlation length exponent is obtained by 
expanding the RG beta functions near the critical fixed 



point. One finds: 

l/u=\ . (40) 

[s + 0{s^) s<l 

These results in excellent agreement with NRG data, as 
shown in Fig. |5| Interestingly, the correlation length ex- 
ponent diverges for both s ^ O"*" and s — > 1~, where the 
second-order transition disappears. Thus the transition 
is bounded by two lower- critical dimensions, somewhat 
similar to the p-h symmetric pseudogap Kondo model de- 
scribed in Sec. I VII However, for s < 1 and r > the two 
problems are in different universality classes^^, as can be 
seen e.g. by comparing the correlation length exponents 
in Figs. and |S1 

For small a, Ao, Eq. (|38|l yields for the phase bound- 
ary 

ac oc AJ"" for Aq <C lOc (41) 

valid for all < s < 1 , which compares well with NRG 
We also quote a few results for the critical exponents 
associated with a local field (|10|l . obtained by the small- 
s expansion of Sec. IVII (]l 

7 = l + 0{s), 
1/S = l-2s + 0{s^). (42) 

Further, one finds the exact result x — y — s which relies 
on the bath bosons being non-interacting^. Thus the 
local spin correlation function C{ljj) diverges as cj^" at 
criticality, whereas it follows uj^ in the disordered phase. 

We note that all critical exponents obey hyperscaling 
which has been verified by NRG; we conclude that the 
critical fixed point is interacting for all < s < 1. The 
exponents (3 and 5 (Fig. I10|l are distinct from the ones 
of the long-range Ising model with l/r^^'^ interaction for 
s < 1/2 which displays mean-field behavior there. The 
obvious failure of the quantum-classical mapping will be 
discussed in Sec. IXIII 

E. Applications 

Spin-boson and dissipative impurity models2i4 have ap- 
plications in many fields like glass physics, mechanic fric- 
tion, damping in electric circuits, electron transfer in bi- 
ological molecules etc. Re-newed interest in spin-boson 
models arises in the quantum computation, for modelling 
the coupling of qubits to a noisy environment and the 
associated decoherence processes. Interestingly, the de- 
scription of 1// noise in electrical cicuits leads to sub- 
ohmic damping with s = (at least over a certain range 
of energies). 

Modifications of standard spin-boson physics include 
the influence of localized modes which interact with two- 
level systems - those modes can be represented by a 
discrete spin system, leading to so-called central spin 
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FIG. 10: Numerical data for the critical exponent 1/5 
{rnj obtained from NRG, together with the RG result (g^l- 
The dashed shows the mean-field value (5 = 3 of the one- 
dimensional Ising model with l/r^'^" long-range interaction. 
After Ref. 0. 



models. Electron transport through quantum dots under 
the influence of external noise or degrees of freedom like 
phonons involves the combined effect of tunneling from 
electronic leads and bosonic damping, with competing 
interaction channels - this naturally leads to Bose-Fermi 
Kondo models described in Sec. IIXI 



SU(2) symmetry is preserved, contrary to the spin-boson 
model (133 • 

In a paramagnetic phase of the bulk, the a, are gapped, 
and the low-energy properties of TYbk are that of a free 
spin. In contrast, at a magnetic bulk critical point in d 
dimensions the Og obey a gapless spectrum of the form 
l)37|l . with exponent s = d — 2. Note that in dimensions 
d < 3, boson self-interactions are relevant in the RG sense 
at the critical point, and then a 0"* theory formulation 
becomes more appropriateiS*ii: 



n 



(b) 



2 4! ^^"'^ 



d^x 



(45) 

where s is the tuning parameter, c a velocite, uq the bulk 
non-linear interaction, and tTc the momentum conjugate 
to 0Q, hence 



ha{x,t),'!Tf3{x',t)] = iSalsS'^ix - x') 



(46) 



VIII. BOSE KONDO MODEL 



In a magnetically ordered bulk phase, the Hamiltonian 
includes the local exchange field from the static order 
parameter, and only transverse magnetic fluctuations re- 
main gapless. 

In the representation H43() of the Bose Kondo problem 
the complexity of the impurity lies in the spin commuta- 
tion relations of Sa the model can also be written using a 
spin coherent state path integral, then the only impurity 
non-linearity appears in the Berry phase termHi^. 



A class of models which is closely related to the 
spin-boson model, termed Bose Kondo models, describe 
magnetic impurities coupled to spinful bosons^. Typ- 
ically these bosons repesent collective bulk spin fluc- 
tuations. The Hamiltonian is conveniently written as 

Hbk -70^a<^a(0). (43) 

In the simplest case, the bosonic bath consists of free 
vector bosons, 

qa 

describing collective spin fluctuations in the d- 
dimensional host material, with a dispersion — + 
(?q^. At zero temperature, the mass to coincides with 
the bulk spin gap A^. A gap Ag > describes a quan- 
tum paramagnet, A^ — corresponds to a bulk quan- 
tum critical point between the paramagnet and an an- 
tiferromagnctically ordered phase; the momentum q is 
measured relative to the ordering wavevector Q. The 
field (/'q(O) in Eq. (|43|) represents the local orientation 
of the antiferromagnetic order parameter, it is given by 
4>aiS^) = Hq^o^qa + ^-qt^l \J ^ql J ■, whcrc J is the bulk ex- 
change constant, i.e., in the parametrization of Eq. H35|) 

— 1/2 

we have Ag cx ojg . Importantly, in this model the 



A. Bose Kondo model: RG near s = 1 

The model (I43|l can be analyzed using a perturbative 
RG expansion in 70 and the bulk interaction uqj which 
are both marginal at s = 1 (d 3). It is important to 
treat both couplings on an equal footing as there is a non- 
trivial "interference" between uq and 70, which appears 
at two-loop order, and the interaction up significantly 
modifies the magnetic environment coupling to the im- 
purity. In the magnetic context it is not permissible to 
treat the environment as a Gaussian quantum noise, and 
focus only on its Kondo-like coupling to the impurity. 

A systematic RG analysis of Wbk was carried ou1ii2*ii, 
resulting in the following RG equation for the renormal- 
ized impurity coupling 7: 

m - -iy^7+7^-7^ + ^ + ^[^(^+l)-l/3] 

(47) 

to two-loop order, and u is the renormalized bulk interac- 
tion. For s < 1, i.e., d < 3, there is a stable intermediate 
coupling fixed point at 7*^ = {1- s)/2 + 0{[l- sf) - this 
"bosonic fluctuating" fixed point (B-FL) displays univer- 
sal local-moment fluctuations as detailed below. Note 
that this is in stark contrast to the model with Ising 
symmetry which flows to strong coupling in the same 
parameter regime - this strong-coupling fixed point is 
apparently prohibited here by the competition between 
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the multiple baths. The B-FL fixed point is unstable 
w.r.t. SU(2) symmetry breaking, but the existence of 
an intermediate-coupling fixed point is preserved for XY 
symmetry. 

For s > 1 the coupling 7 is (marginally) irrelevant. 
Here, the low-tcmpcrature properties of TYbk are equiv- 
alent to an asymptotically free spin. 

B. Bose Kondo model: RG near s = —1 

An alternative approach to the dynamics of an impu- 
rity embedded in a bulk magnet is provided by a repre- 
sentation of the bulk spin fluctuations using a non-linear 
sigma model. The fluctuations of the order parameter 
amplitude, become irrelevant in low dimensions, and 
a hard-spin representation by a unit-length held Na {x,t), 
capturing only angular fluctuations of the Neel order pa- 
rameter, is appropriate. The non-linear sigma model for- 
mulation allows for an RG expansion of the bulk quantum 
critical properties in (l-f-e) dimensions. 

The Bose Kondo problem can be tackled using such a 
fixed-length representation^. Interestingly, in the scal- 
ing limit the quantum impurity behaves as if it is in the 
7 — > 00 limit, and hence the universal nature of the cou- 
pling between the bulk and impurity is explicit. The 
properties of the stable impurity fixed point at bulk crit- 
icality can be computed in a systematic expansion in 
{d — 1) (i.e., near s — —1 in terms of the bath spectrum 
exponent). All results were found to be consistent with 
those of the (3 — d) expansion of the previous subsection. 

C. Bose Kondo model: large- A'^ 

Generalizing the impurity spin to SU(iV) symmetry 
allows for a large- limit corresponding to a dynamic 
saddle point With the impurity representation 125|) 
the Hamiltonian for the impurity takes the form 

-HbK = E \fiMb,aP + + Ao/l/a ■ (48) 

For N ^ oo we obtain an integral equation for the / 
fermion Green's function: 

S/(r) = 4'^)(r)G/(r) (49) 

where Gq'^'' is the local Green's function of the bath, 
G^o\r) = -E„A2(T.(&, + fet_^)(r)(&t +6_^)(0)), and 
the self-energy S/ is defined by: 

"^(iw„) = iw„ - Aq - I]/(ia;„) . (50) 

These equations correspond to the summation of all 
self-energy diagrams with non-crossing bath boson lines, 
and can be understood as the bosonic analogue of the 
fermionic NCA equations (|30I31I32|I . 



The equations (|49I50|I have been analyzed in Ref. ITlI 
One obtains a G/ solution with a power-law behavior at 
low energies, corresponding to an intermediate-coupling 
fixed point. Many of the properties are similar to the 
SU(2) case with non-interacting bosons, analyzed above 
using RG. Interestingly, the NCA-type approach can be 
extended to a finite concentration of impurities and al- 
lows to capture the feedback of the impurity scattering 
to the bulk spin excitationsii. 



D. Results for observables 

Although the Bose Kondo model does not display a 
quantum phase transition, we quote a few results for 
the intermediate-coupling fixed point obtained from the 
RG expansion. Some of the properties for interacting 
bosons in d = 2 have been also investigated by quan- 
tum Monte Carlo simulations of two-dimensional lattice 
antiferromagnet3^2iffi. 

The structure of the fixed point implies that the impu- 
rity shows universal local- moment fluctuations: the lo- 
cal susceptibility xioc obeys a scaling form H12|l with an 
anomalous exponent 77^^ > for s < 1, which means that 
local spin correlations at T = are characterized by a 
power lawifliSiiSS, {S{t)S) oc r^''^. For non-interacting 
bath bosons, one obtains the exact result 77^ — 1 — s 
[or X — y = s according to ifTUIl ]. Interestingly, this ex- 
act exponent implies that the impurity fluctuates faster 
than its environment for s < 0, which violates a theorem 
due to Griffiths, stating that the impurity has to fluctu- 
ate slower than the environment, 77^ < 1 + s. Thus, the 
Bose Kondo problem with non-interacting bosons likely 
changes qualitatively at s = 0, however, details are un- 
known. In contrast, higher powers of (1 — s) appear in 
the perturbative expansion of 77^ for interacting bosons, 
and the results are expected to fulflU the above inequality 
for aU -1 < s < 1. 

Other observables at the intermediate-coupling flxed 
point can be associated with the properties of a fractional 
spin: The Curie contribution to the impurity susceptibil- 
ity is Tximp = Cimp, and the impurity entropy S'imp is a 
finite constant and larger than In 2. Both Cimp and S'imp 
are universal functions of s (or the dimensionality d of 
the bath), but depend on whether the bosons are inter- 
acting and non-interacting. Related theoretical results 
were obtained recently^* in magnetically ordered states 
in the presence of spin anisotropy, and in Ref. Is^ . 

We summarize that the properties of the Bose Kondo 
problem with a bath given by a quantum critical magnet, 
i.e., with interacting bosons, appear to vary smoothly for 
1 < c? < 3 (i.e., — 1 < s < 1). For d — 1 the problem is 
ill-defined as the bulk transition disappears (and the en- 
vironment at the putative transition ceases to fluctuate), 
and for d > 3 the impurity coupling is (marginally) irrele- 
vant, resulting in an asymptotically free spin. In contrast, 
in the Bose Kondo model with a non-interacting (Gaus- 
sian) bath the impurity behavior is believed to change 
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qualitatively at s = 0, i.e., non-trivial intermediate cou- 
pling physics is only present for < s < 1. 

E. Applications 

Vector baths naturally appear in the context of bulk 
spin fluctuations. Thus impurity moments embedded in 
magnets are naturally described by Bose Kondo models - 
such models become appropriate when the environment 
of the impurity is in the vicinity of a magnetic ordering 
transition, and there are low-energy spin excitations in 
the bulk; the latter may be viewed as excitonic particle- 
hole bound states of a metal/insulator/superconductor 
which peel off below the continuum of a p air of f ermionic 
particles or holes. As discussed in Refs. llOlllI the self- 
interaction of the bulk bosons cannot be neglected be- 
low the upper-critical dimension of the bulk, as it also 
changes the impurity physics. Concrete realizations are 
magnetic atoms in insulating bulk magnets. At low en- 
ergies, these considerations can be extended to d-wave 
superconductors, and can describe non-magnetic impu- 
rities like Zn or Li in cuprates, which are known to in- 
duce magnetic moments. For a finite impurity concen- 
tration, the universal interaction between impurity mo- 
ments and host spin fluctuations leads to universal im- 
purity damping of spin fluctuations in cuprate supercon- 
ductors - corresponding signatures have been observed in 
inelastic neutron scattering in Zn-doped YBa2Cu307^^. 

Ref. 1571 discussed the photoemission spectrum of a con- 
duction electron in an insulator in the vicinity of the mag- 
netic transition^. Let us focus on the minimum energy 
of its dispersion. While away from the critical point a 
sharp quasiparticle pole appears, the interaction with the 
magnetic modes leads to a universal damping of this pole 
at the magnetic bulk transition, resulting in a power-law 
spectrum G{uj) ^ l/(w — £o)^ - This behavior reflects 
an orthogonality catastrophy in the spin sector, and is 
the bosonic analogue of the familiar X-ray edge problem 
in a Fermi liquid^^. 

Spins coupled to bulk spin fluctuations can also occur 
in quantum dots coupled to leads with magnetic collec- 
tive modes; taking into account the fermionic degrees of 
freedom as well leads to Bose-Fermi Kondo models as 
described in Sec. IIXI 

IX. BOSE-FERMI KONDO MODEL 

Novel phenomena occur for magnetic impurities cou- 
pled to both a fermionic and a bosonic bath, where the 
bosons represent e.g. collective host spin excitations. In 
the resulting Bose-Fermi Kondo model the two interac- 
tions compete in a non-trivial manneriSiiS^, and fermionic 
Kondo screening can be strongly suppressed by host spin 
fluctuations. The Hamiltonian is a straightforward com- 
bination of the bosonic and fermionic Kondo models, 
Ti. = Hbfk + ^b''' + with an impurity spin 5 = 1/2 



and 

WbFK = JaSa<Pa{0) + JkS'qSq(O) (51) 

with notations as above. The spin-1 order parameter field 
(j)a and the spin-1/2 fermions Cka, with local spin density 

Sa(0), are described by the baths H^^^ and 'H'^\ respec- 
tively. Most interesting is again the case of a bosonic bath 
with zero or small gap, corresponding to the vicinity to 
a magnetic quantum critical point in the bulk. 



A. Bose-Fermi Kondo model: Phase diagram and 

RG 

We now discuss the phase diagram of the Bose-Fermi 
Kondo model, based on perturbative RG and large-TV 
results. We restrict ourselves here to a metallic density of 
states for the fermions (the pseudogap case was discussed 
in Ref.El. 

For a gapless bosonic bath with spectrum uj", s < 1, 
corresponding to a magnetic critical point in 2 -f s di- 
mensions, both the fermionic and bosonic couplings are 
relevant, and compete. For large Jk fermionic Kondo 
screening wins, resulting in a fully screened spin. Large 
7o can completely suppress Kondo screening, driving the 
system into the intermediate-coupling fixed point of the 
Bose-Kondo problem, with universal local-moment fluc- 
tuations (Sec. IVIIII The competition is captured by the 
RG equationaiSia^: 

m = -f+n^- (52) 

The phase diagram for the Bose-Fermi Kondo model 
thus shows a Kondo-screened phase, a bosonic fluctu- 
ating phase, and a continuous quantum phase transition 
in between. The RG flow is shown in Fig. ^] Both 
intermediate-coupling fixed points are perturbatively ac- 
cessible for small (1 — s); it is likely that the structure 
of the phase diagram also applies to s = (i.e. d = 2), 
however, no accurate numerical calculations are available 
to date. (As with the Bose Kondo model, for interact- 
ing bath bosons the properties likely change smoothly for 
— 1 < s < 1, whereas in the Gaussian case a qualitative 
change is expected for s < 0.) Both the critical and the 
bosonic fluctuating fixed points are unstable w.r.t. break- 
ing of SU(2) symmetry, but the structure of the phase di- 
agram is similar for both XY and Ising symmetries, with 
the difference that in the Ising case the boson-dominated 
phase corresponds to a strong-coupling fixed pointS. 

If the bosonic bath is gapped, or its spectrum has 
exponent s > 1, then the coupling to the bosons is 
(marginally) irrelevant, and the low-temperature prop- 
erties of (|51ll are that of the familiar fermionic Kondo 
model. However, the Kondo temperature can be strongly 
renormalized due to the competition from the bosonic 
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FIG. 11: RG flow diagram for the Bose- Fermi Kendo 
model. The two axes denote the renormahzed couphngs to the 
fermionic (j) and bosonic (7) baths, respectively. LM is the 
local-moment fixed point of a decoupled impurity, F-SC the 
fermionic Kondo-screened fixed point, B-FL the intermediate- 
coupling fixed point described in Sec. lVllH and FB-CR is crit- 
ical fixed point of the Bose-Fermi Kondo problem. A metallic 
density of states of the fermionic bath is assumed, and the 
bosonic bath has s < 1, i.e., d < 3 in the magnetic context. 
Aft^er Ref. HI 



bath. An exception is the case of s = 1 and Ising sym- 
metry: here the bosonic coupling drives a Kosterlitz- 
Thouless transition between a screened and unscreened 
impurity^^. 



B. Bose-Fermi Kondo model: large- 

A multi-channel version of the Bose-Fermi Kondo 
model can be analyzed in a large-A'' limit, with a an- 
tisymmetric representation of SU(A'^) for the impurity 
spir^^. The bosonic bath is treated as described for 
the Bose Kondo model fSec. IVIII C|) . and the fermionic 
bath is captured by the standard multi-channel NCA 
approach^^. One obtains a set of coupled integral equa- 
tions where the self-energy of the auxiliary fermions now 
has contributions from both the bosonic and fermionic 
baths. In the low-energy limit power-law solutions for 
the auxiliary-particle propagators appear, corresponding 
to intermediate-coupling physics. 

The phase diagram2& is similar to the one obtained 
from the perturbative RG described above, with the 
difference that large Kondo coupling produces a stable 
intermediate-coupling fixed point corresponding to over- 
screening (in contrast to the exactly screened situation 
of the single-channel model). Importantly, oj/T scal- 
ing (|12|l is obtained for all values of the bath exponent 
-1< s < 1. 



C. Results for observables 

The properties of the stable phases of the Bose- 
Fermi Kondo model have already been discussed, see 
Sec. IVIII Dl for the boson-dominated B-FL fixed point. 
Notably, magnetic properties of the quantum critical 
point of (|51|l are very to the B-FL fixed poinljSS: the 
local susceptibility follows the scaling law H12|l with an 
anomalous exponent 7/^ > for s < 1, and the impu- 
rity susceptibility features a fractional Curie moment. In 
addition, the conduction electron T matrix at criticality, 
discussed in Ref. is (only) logarithmically singular for 
a metallic fermionic bath. We note, however, that reli- 
able numerical results for low temperatures are lacking. 

D. Applications 

In the context of strongly correlated electron sys- 
tems, which often feature Fermi-liquid quasiparticles and 
strong spin fluctuations at the same time, the question 
of the interplay between fermionic and bosonic Kondo 
physics, as described by the Bose-Fermi Kondo model, 
arises. It has been proposed^ that this interplay plays 
a role in high- Tc cuprates, where NMR experiments in- 
dicate that dilute impurity moments are screened at op- 
timal doping, but Tk is essentially suppressed to zero in 
underdoped compounds^''. A Fermi-Bose Kondo model, 
taking into account both the pseudogap density of states 
of the Bogoliubov quasiparticles and the strong anti- 
ferromagnetic fluctuations, provides a natural explana- 
tion: spin fluctuations increase with underdoping, thus 
strongly reducing Tk due to the vicinity to the boundary 
transition which exists even in the absence of a bosonic 
bathes. 

We note that other approaches to Kondo physics in 
correlated environments have been put forward, which 
do not necessarily rely on the heuristic separation of 
bosonic and fermionic bulk degrees of freedom. Near 
a magnetic quantum critical point, ordered magnetiza- 
tion droplets form around impurities and interact with 
the electronic environment^^. A related work on clus- 
ter Kondo physics^- used an explicit microscopic model 
for a magnetic nanostructure coupled to a metallic bath. 
Further, it has been pointed out^ that long-range criti- 
cal bulk spin correlations can open additional screening 
channels (even for an initially local impurity), resulting 
in multi-channel Kondo physics at intermediate tempera- 
tures. The class of cluster Kondo problems in disordered 
correlated environments is closely related to the issue of 
quantum Griffiths behavior near magnetic quantum crit- 
ical points, with possible realizations in chemically sub- 
stituted heavy-fermion metals. 

Various modifications of Bose-Fermi Kondo physics ap- 
peared recently in the context of mesoscopic systems. 
Quantum dots coupled to ferromagnetic leads involve 
both fermionic and bosonic baths plus an additional mag- 
netic field resulting arising from the lead polarization. 
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Charge fluctuations on a metallic island, described by 
Matveev's charge Kondo effect'^'^, can be significantly 
modified by a noisy electromagnetic environment, with 
the low-energy physics described by an anisotropic Bosc- 
Fcrmi Kondo model^'^. 

Finally, we mention that the Bose-Fermi Kondo model 
has recently received a lot of interest in the con- 
text of extended dynamical mean-field theory^® where 
a lattice model is mapped onto a self-consistent im- 
purity model with both fermionic and bosonic baths. 
Motivated by neutron scattering experimental^ on the 
heavy-fermion compound CeCug-xAux which indicate 
momentum-independent critical dynamics at an antifer- 
romagnetic ordering transition, a self-consistent version 
of the Bose-Fermi Kondo model has been proposed to 
describe such "local" critical behavior. In this scenario, 
the critical point of the lattice model is mapped onto the 
FB-CR critical fixed point of the impurity model. 

X. TWO-IMPURITY KONDO MODELS 

Models of multiple impurities offer a new ingredi- 
ent, namely the exchange interaction between the im- 
purity spins, which competes with Kondo screening of 
the individual impurities. This inter-impurity interac- 
tion, which can lead to a magnetic ordering transition 
in lattice models, arises both from direct exchange and 
from the Ruderman-Kittel-Kasuya-Yosida (RKKY) in- 
teraction mediated by the conduction electrons. 

Let us consider a Kondo model with two impurity spins 
1/2 Si, S2, and a Hamiltonian written as 

^ = II^S+II"^kS, •s,(0) + Hi2. (53) 

i i 

Each impurity spin is coupled to its individual bath Tib.i, 
i.e., there is a total of two screening channels, and 'H12 
contains the inter-impurity interaction. Such models of 
two coupled impurities or two-level systems have been 
investigated in a number of papers; the resulting phases 
and phase transitions depend on the symmetry properties 
of the bath and the couplings, and will be discussed in 
the following. 

A. Behavior for SU(2)-symmetric coupling 

The case of SU(2)-symmetric direct exchange coupling 
between two impurity S = ^ spins, 

Hi2 = ifSi-S2, (54) 

has been investigated in detail in the literatureSSiiiii^. 
Two different regimes are possible as function of the 
inter-impurity exchange K: for large antiferromagnetic 
K the impurities combine to a singlet, and the interac- 
tion with the conduction band is weak, whereas for fer- 
romagnetic K the impurity spins add up and are Kondo- 
screened by conduction electrons in the low-temperature 



limit. It has been shown that these two parameter 
regimes are continuously connected (without a T = 
phase transition) as K is varied in the generic situation 
without particle-hole symmetry. Notably, in the particle- 
hole symmetric case one finds a transition associated with 
an unstable non- Fermi liquid fixed poin"fe2ii2^. This fixed 
point is somewhat similar to the one of the two-channel 
Kondo model of Sec. EI e.g., it features a residual entropy 
of i In 2. 

We note that different physics is encountered by cou- 
pling the two-impurity system to a single conduction 
band channel onljsS. For ferromagnetic K complete 
Kondo screening by the conduction electrons is prohib- 
ited, resulting in an underscreened Kondo phase. Upon 
variation oi K a, Kosterlitz-Thouless-type transition be- 
tween a singlet and a doublet state occurs, associated 
with a second exponentially small energy scale in the 
Kondo regime^^. 

The physics becomes even richer if multi-channel 
physics is combined with multi-impurity physics - here, 
a variety of fixed points including such with local non- 
Fermi liquid behavior can be realized. 

B. Behavior for Ising coupling 

The two- impurity model with Ising inter-impurity cou- 
pling, 

H12 = KJlS'^ , (55) 

is qualitatively different from the SU(2)-symmetric one^^. 
It is particularly interesting, because for large \Kz\ the 
two Ising-coupled spins form a magnetic mini-domain 
which still contains an internal degree of freedom as the 
ground state of H12 is doubly degenerate (in contrast 
to the inter-impurity singlet mentioned above). For the 
case of antiferromagnetic the two low-energy states 
of the impurities (forming a pseudospin) are | Ti) ^J^d 
I If). The fate of this pseudospin degree of freedom 
depends on the strength and asymmetry of the Kondo 
coupling J between the spins and the bath electrons 
On the one hand, for small Kondo couplings Jz and 
large if^, tunneling between the two pseudospin config- 
urations, I tJ,) and I |t)i is supressed at low energies, 
i.e., the mini-domain is "frozen" as T — > 0, and the 
ground state entropy is 5*0 = In 2. On the other hand, 
for small the two impurities are individually Kondo 
screened, resulting in a Fermi- liquid phase with vanishing 
residual entropy. A quantum phase transition occurs for 
Kz ^ where T^p is the single-impurity Kondo tem- 
perature. This phase transition is of Kosterlitz-Thouless 
type and can also be tuned by varying the Ising compo- 
nent Jz of the Kondo coupling. This leads to the pres- 
ence of a second exponentially small scale T*, which is 
a collective Kondo temperature associated to pseudospin 
screening. The high-temperature In 4 impurity entropy is 
quenched in two stages: first, at the scale T° « Kz^ the 
mini-domain "forms" , quenching half of entropy; second 
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the strong fluctuations kill the remaining In 2 entropy at 
the much lower scale T*— . 



C. Applications 

Experimentally, quantum dots provide an ideal labo- 
ratory to study systems of two (or more) "impurities" 
- note that the local "impurity" states can arise either 
from spin or charge degrees of freedom on each quantum 
dot. 

In particular, a number of experiments have been 
performed on multi-level or coupled quantum dot sys- 
tems which can be directly mapped onto models of two 
Kondo or Anderson impurities^^. Experiments reported 
in Ref. have used a setup of two small quantum dots 
coupled to a larger metallic island and various external 
leads. Varying gate voltages allowed to tune the spin 
states of the dots as well as the RKKY inter-dot coupling, 
and signatures of the transition discussed in Sec. IX Al 
should be observable. 

If charge degrees of freedom are employed for the 
Kondo effect^ then two capacitively coupled dots present 
a promising realization of the Ising coupling of Sec. IX Bl 
A Kosterlitz-Thouless quantum phase transition will oc- 
cur upon variation of the capacitive inter-dot coupling. 
Adding a small electron tunneling allows transport exper- 
iments through this setup, and a universal conductance 
anomaly near the transition has been predicted2&. 

Using recent advances in scanning probe techniques, 
experimental realizations of multi-impurity models us- 
ing magnetic adatoms on metallic surfaces are possibleiS. 
Here, the spectral functions of the inidvidual impurities 
can be spectroscopically measured, and the suppression 
of the Kondo effect for impurity dimers has been demon- 
strated. One interesting set of results is concerned with 
three Cr adatoms forming a trimer - this system has 
been proposed to show a stable fixed point with local non- 
Fermi liquid behavior, provided that the threefold spatial 
symmetry is unbrokenSi. Cluster of magnetic impurities, 
i.e., magnetic nanostructures, embedded in a metallic en- 
vironment represent a natural extension of these lines of 
thought. An initial study of such a Kondo "necklace" 
model^^ showed that multi-channel Kondo physics can 
be realized in a sizeable temperature window. 

XI. MULTI-ORBITAL ANDERSON MODELS 

Similar to coupled impurities, multi-orbital or multi- 
level impurities offer a larger set of local degrees of free- 
dom which can undergo a boundary phase transition. 

Experimental realizations can be impurity atoms with 
multiple low- lying crystal field states, offering bath 
charge (e.g. quadrupolar) and spin degrees of freedom. 
In the presence of orbital and spin degeneracy the result- 
ing model can show e.g. quadrupolar two-channel Kondo 
effect (see Sec.0). A distinct but related non-trivial fixed 



point has been found for a two-orbital Anderson model 
with Hund's rule coupling^. 

Alternatively, multi-level quantum dots can be tuned 
to study Kondo physics beyond the simple Kondo model. 
For instance, for two impurity levels and two electrons 
in the low-energy sector the impurity can be close to a 
singlet-triplet transition, which leads to a strongly en- 
hanced Kondo temperature. Other novel phenomena like 
a SU(4) Kondo effect have also been discussed. One can 
envision various quantum phase transitions in those de- 
vices which have not yet been fully explored. 

XII. QUANTUM-CLASSICAL MAPPING 

Significant progress in impurity physics has been made 
on the basis of various mappings. On the one hand, 
suitable transformations (e.g. bosonization, flow equa- 
tions) can lead to great simplifications. On the other 
hand, Gaussian (non-interacting) bath degrees of free- 
dom can be formally integrated out, leading to a long- 
range self-interaction in imaginary time direction for the 
impurity. Discretizing imaginary time then leads to 
one-dimensional statistical mechanics models with long- 
range interactions - this is a special case of the standard 
quantum-classical mapping for quantum phase transi- 
tions. (For bath particles with a relevant self-interaction 
the idea of the quantum-classical mapping is inapplica- 
ble in general, as the bath cannot be formally integrated 
out.) 

A. Ising symmetry 

Specifically, the spin-boson model fSec lVIlT) with bath 
density of states yields an effective impurity self- 
interaction 

5i„t = J dTdT'a,{T)g{T - r')a,(T') (56) 

with ^(t) cx l/r^+* at long times. The classical counter- 
part is the one-dimensional Ising mode&^ 

Ki = -^J»,5f5| + HsR (57) 

(v) 

with interaction — J/\i — TisR contains an ad- 

ditional generic short-range interaction which arises from 
the transverse field, but is believed to be irrelevant for 
the critical behavior*^. 

For s = 1 the quantum-classical mapping appears 
consistent: both the ohmic spin-boson model and the 
1/r^ Ising model show a Kosterlitz-Thouless transi- 
tion. Notably, also the metallic Kondo model falls in 
this universality class, which can be established using 
bosonizationpSi^. 

The three models (classical Ising, anisotropic Kondo, 
spin-boson) display non-trivial second-order transitions if 
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the characteristic bath exponent is varied. However, the 
transitions of the fermionic and bosonic quantum models 
are no longer in the same universality class! This has 
been shown by explicit calculations, see Sec. IVII for the 
pseudogap Kondo model and Sec. IVIII for the sub-ohmic 
spin-boson model. For instance, the particle-hole pseu- 
dogap Kondo model shows a phase transition only for 
< 7' < i (and has a complicated fixed point structure 
in the asymmetric case, with a transition for all r > 0), 
whereas the spin-boson model has a phase transition for 
< s < 1. These fundamental differences arise from the 
distinct character of the bath degrees of freedom: Lin- 
early dispersing fermions in (1 -I- r) dimensions cannot be 
bosonized, and a direct evaluation of the corresponding 
fermionic determinants is not easily possible. Contrary 
to metallic Kondo problems, the pseudogap Kondo model 
is not conformally invariant, and conformal field theory 
techniques cannot be applied. 

The bosonic case appears simpler, and the formal map- 
ping between the spin-boson model and the classical 
Ising model H57|) is straightforward. Indeed, the spin- 
boson model and the classical Ising model remain equiv- 
alent in the range ^ < s < 1. However, in Ref. |9l it 
has been shown that this equivalence breaks down for 
< s < i; The classical model is effectively above its 
upper-critical dimension--, given hy d = 2s, the tran- 
sition takes mean-field exponents, and hyperscaling is 
violated. The spin-boson model, however, displays an 
interacting fixed point with non-trivial exponents fulfill- 
ing hyperscaling, see Sec. IVIII How can the quantum- 
classical mapping fail? To establish the mapping a Trot- 
ter decomposition of the quantum partition function is 
employed where the imaginary axis of length /3 = 1/T 
is divided into N slices, leading to an Ising chain H57() 
with N sites. This procedure is exact when the limits 
N oo and (3 —^ oo are taken in this order. However, 
the limit f3/N — > leads to a diverging near-neighbor cou- 
pling in the term Hsr of the classical Ising model (|57pi. 
This may in fact change the critical behavior of Hd (|57|l . 
In other words, the quantum and classical problems are 
only equivalent if the low-energy behavior of Tid is inde- 
pendent of the order of limitg^i. Apparently, these two 
limits cannot be interchanged for s < 1/2, and the naive 
quantum-classical mapping fails to describe the critical 
properties of the strongly sub-ohmic spin-boson model. 



B. XY and SU(2) symmetries 

For bosonic quantum impurity models with continuous 
symmetry our main knowledge arises from the RG cal- 
culations sketched in Sec. lVIlTI For a quantum spin cou- 
pled to non-interacting vector bosons the naive quantum- 
classical mapping would predict the equivalence to a one- 
dimensional classical XY or Heisenberg model with long- 
range interactions. These models are knows to display 
conventional ordered and disordered phases. 

In contrast, the RG analysis of Sec. I Villi uncovered the 



existence of a stable intermediate-coupling fixed point 
with non-trivial power laws; the inclusion of an external 
field will drive a transition to a disordered phase. Thus, 
for XY and SU(2) symmetries even the stable (ordered) 
phase of the quantum model is not correctly reproduced 
by the quantum-classical mapping for any value of s < 1, 
and the same likely applies to the transition properties. 
A related observation appears in a SU(iV) version of the 
Bose-Fermi Kondo modetf^ where uj/T scaling is obeyed 
for all s < 1, contrary to the expectations from the clas- 
sical model. 

It is tempting to associate the general inapplicability 
of the mapping with the presence of the Berry phase term 
describing the impurity spin dynamics - this Berry phase 
is imaginary in the path integral representation and has 
no classical analogue. Thus we conjecture that the map- 
ping works for a quantum rotor coupled to a Gaussian 
bath, where the Berry phase term is absent. 



XIII. SUMMARY 

We have reviewed a variety of zero-temperature critical 
points in quantum impurity models, both with fermionic 
and bosonic baths. Significant progress has been made 
in recent years, both in analytical and numerical work, 
which has uncovered e.g. the existence of both lower- 
critical and upper-critical dimensions and associated per- 
turbative epsilon-expansions. In addition, the validity of 
the naive quantum-classical mapping between quantum 
impurity models and one-dimensional statistical mechan- 
ics models with long-range interactions has been critically 
examined. 

The impurity quantum phase transition have a variety 
of applications in single-impurity models, e.g., for im- 
purities embedded in correlated electronic environments, 
and in quantum dot systems with potential applications 
for spintronics. In addition, effective impurity models 
appear in dynamical mean-field approximations to corre- 
lated lattice models where they are supplemented by self- 
consistency conditionsfii. Intermediate-coupling fixed 
points of impurity models have been proposed to describe 
unconventional "liquid" phases of spin glass models^ and 
so-called local quantum critical points®*. 

Interesting future directions include the study of non- 
equilibrium properties near equilibrium quantum phase 
transitions, which are relevant, e.g., for quantum dot 
transport experiments and for quantum decoherence pro- 
cesses. 
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